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Abstract. We explore the recently introduced statefinder parameters. After reviewing their basic properties, 
we calculate the statefinder parameters for a variety of cosmological models, and investigate their usefulness as 
a means of theoretical classification of dark energy models. We then go on to consider their use in obtaining 
constraints on dark energy from present and future supernovae type la data sets. Provided that the statefinder 
parameters can be extracted unambiguously from the data, they give a good visual impression of where the correct 
model should lie. However, it is non-trivial to extract the statefinders from the data in a model-independent way, 
and one of our results indicates that an expansion of the dark energy density to second order in the redshift 
is inadequate. Hence, while a useful theoretical and visual tool, applying the statefinders to observations is not 
straightforward. 
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1. Introduction 

It is generally accepted that we live in an accelerat- 
ing universe. Early indications of this fact came from 
the magnitude-redshift relationship of galaxies (Solheim 
1966), but the reality of cosmic acceleration was not taken 
seriously until the magnitude-redshift relationship was 
measured recently using high-redshift supernovae type la 
(SNIa) as standard candles (Riess et al. 1998, Perlmutter 
et al. 1999). The observations can be explained by in- 
voking a contribution to the energy density with negative 
pressure, the simplest possibility being Lorentz Invariant 
Vacuum Energy (LIVE), represented by a cosmological 
constant. Independent evidence for a non-standard con- 
tribution to the energy budget of the universe comes from 
e.g. the combination of the power spectrum of the cosmic 
microwave background (CMB) temperature anisotropies 
and large-scale structure: the position of the first peak 
in the CMB power spectrum is consistent with the uni- 
verse having zero spatial curvature, which means that the 
energy density is equal to the critical density. However, 
several probes of the large-scale matter distribution show 
that the contribution of standard sources of energy den- 
sity, whether luminous or dark, is only a fraction of the 
critical density. Thus, an extra, unknown component is 
needed to explain the observations (Efstathiou et al. 2002; 
Tegmarket al. 2004). 



Several models describing an accelerated universe have 
been suggested. Typically, they are tested against the 
SNIa data on a model-by-model basis using the rela- 
tionship between luminosity distance and redshift, di^(z), 
defined by the model. Another popular approach is to 
parametrize classes of dark energy models by their pre- 
diction for the so-called equation of state w{z) = Px/ p^, 
where px and px are the pressure and the energy den- 
sity, respectively, of the dark energy component in the 
model. One can then Taylor expand w{z) around z = 0. 
The current data allow only relatively weak constraints 
on the zeroth-order term wq to be derived. A prob- 
lem with this approach is that some attempts at ex- 
plaining the accelerating Universe do not involve a dark 
component at all, but rather propose modifications of 
the Friedmann equations (Deffayet 2001; Deffayet, Dvali 
& Gabadadze 2002; DvaU, Gabadadze & Porrati 2000; 
Freese & Lewis 2002; Gondolo & Freese 2003; Sahni 
& Shtanov 2003). Furthermore, it is possible for two 
different dark energy models to give the same equa- 
tion of state, as discussed by [Padmanabhan (2002) and 
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Padmanabhan fc Choudhury (2003']] " 

Recently, an alternative way of classifying dark energy 
models using geometrical quantities was proposed (Sahni 
et al. 2003, Alam et al. 2003). These so-called statefinder 
parameters are constructed from the Hubble parameter 
H{z) and its derivatives, and in order to extract these 
quantities in a model-independent way from the data, one 
has to parametrize H in an appropriate way. This ap- 
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proach was investigated at length in Alam et al. (2003) us- 
ing simulated data from a SNAP^-type experiment. In this 
paper, we present a further investigation of this formal- 
ism. We generalize the formalism to universe models with 
spatial curvature in Section 2, and give expressions for 
the statefinder parameters in several specific dark energy 
models. In the same section, we also take a detailed look 
at how the statefinder parameters behave for quintessence 
models, and show that some of the statements about 
these models in [Alam et al. (2003) have to be modified. 
In Section 3 we discuss what can be learned from current 
SNIa data, considering both direct y^ fitting of model pa- 
rameters to data, and statefinder parameters. In Section 4 
we look at simulated data from an idealized SNIa survey, 
showing that reconstruction of the statefinder parameters 
from data is likely to be non-trivial. Finally, Section 5 
contains our conclusions. 

2. Statefinder parameters: definitions and 
properties 

The Friedmann-Robertson- Walker models of the universe 
have earlier been characterized by the Hubble parameter 
and the deceleration parameter, which depend on the first 
and second derivatives of the scale factor, respectively: 
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where dots denote differentiation with respect to time t. 
The proposed SNAP satellite will provide accurate deter- 
minations of the luminosity distance and redshift of more 
than 2000 supernovae of type la. These data will permit a 
very precise determination of a{z). It will then be impor- 
tant to include also the third derivative of the scale factor 
in our characterization of different universe models. 

Sahni and coworkers (Sahni et al. 2003, Alam et al. 
2003) recently proposed a new pair of parameters (r, s) 
called statefinders as a means of distinguishing between 
different types of dark energy. The statefinders were in- 
troduced to characterize flat universe models with cold 
matter (dust) and dark energy. They were defined as 



(3) 
(4) 



Introducing the cosmic redshift 1 + z = l/a = x, we have 
H = —H'H/a, where H' = dH/dx, the deceleration pa- 
rameter is given by 

qix)^—x-l. (5) 

Calculating r, making use of a' = —a? , we obtain 

rix)^l-2-x+ij^ + — ]x^ (6) 
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^ see _http://snap. lbl.gov 



The statefinder s(a;), for flat universe models, is then 
found by inserting the expressions lO and Q into equa- 
tion Q). The generalization to non-flat models will be 
given in the next subsection. 

The Friedmann equation takes the form^ 



H' 



8t:G 



{Pn 



Pk) 



(7) 



3 a^ 

where pm is the density of cold matter and px is the den- 
sity of the dark energy, and k = —1,0, 1 is the curvature 
parameter with fc = corresponding to a spatially flat uni- 
verse. The dust component is pressureless, so the equation 
of energy conservation implies 

(8) 



Pra = 


PmOC 


I . 






This 


gives 


for the density of dark energy: 


Px = 


Pc- 


Pm 


3fc 
87rGa2 




= 


3 


-AH' 


+ fca;2-anoi?oa^'), 





(9) 



where and fimo and fixo are the present densities of mat- 
ter and dark energy, respectively, in units of the present 
critical density pco — 3Hq/8ttG. In the following, we will 
use the notation Hi = 8TTGpi{t)/H'^{t), flio = ili(i = to), 
where to is the present age of the Universe, and also 
ri = X^i ^i- From Friedmann's acceleration equation 
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where pi is the contribution to the pressure from compo- 
nent i, it follows that 
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Hence, if dark energy is described by an equation of state 
Px — w{x)px, we have 
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In the following subsections, we calculate statefinder 
parameters for universe models with different types of 
dark energy. 

2.1. Models with an equation of state p = w{z)p 

First we consider dark energy obeying an equation of state 
of the form p^ = wp^, where w may be time-dependent. 
Quintessence models (Wetterich 1988; Peebles & Ratra 
1988), where the dark energy is provided by a scalar field 
evolving in time, fall in this category. The formalism in 
Sahni et al. (2003)| and [Alam et al. (2003)| will be gener- 
alized to permit universe models with spatial curvature. 
Then equation Q is generalized to 

r — Q, 

'=3(g-r>/2)' ^^^^ 



^ Throughout this paper we use units where the speed of 
hght c = 1. 
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where fl — flm + ^^ — 1 — ^k, and fik — —k/{a'^H^). 
The deceleration parameter can be expressed as 



1 



1 



[a„ + (1 + 3w)Q^] = -{n + 3wn^) 



(14) 



Differentiation of equation Q together with equation ^ 
leads to 
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From equation (|14|l we have 
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and 

ij = -i/2(l + q), 

giving 

Pc = -2(l + q)i?pc, 
which leads to 

n^ — + 2(l + q)Hn. 
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For cold matter, p,n = —SHpm, giving 

i7,„ = (2g - l)i/r!„„ (22) 

and for the dark energy, p^ = —3(1 + w)Hp^, giving 

n^ = {2q-l-Sw)Hn^. (23) 

Inserting equations H22|) and H23|) into equation H16() and 
the resulting expression into (|15|l finally leads to 



O. 



l + -w{l + w) 



n^~^j^n^. (24) 

Inserting equation H24|) into equation p3|l gives 
s = l + w-^— ^. (25) 

fix = 1 and the expression for 



1 w 
~ 3^' 

For a flat universe H. 
r simplifies to 



9 /-, x^ 3 w^ 

r = 1 + -w(l + w)Q.j^ - -— f^x- 



Note that for the case of LIVE, w = —1 = constant, and 
one finds r = fJ, s = for all redshifts. For a model with 
curvature and matter only one gets r = 2q = f^m, s — 2/3. 
The same result is obtained for a flat model with matter 
and dark energy with a constant equation of state w = 



— 1/3, which is the equation of state of a frustrated net- 
work of non-Abelian cosmic strings (Eichler 1996; Bucher 
& Spergel 1999). Thus, the statefinder parameters cannot 
distinguish between these two models. However, neither 
of these two model universes are favoured by the current 
data (for one thing, they are both decelerating), so this is 
probably an example of academic interest only. 

For a constant w, and llmo + ^xo = li the q-r plane 
for different values of f^x and w is shown in figure ^ 
Quintessence with w = constant is called quiessence. The 
relation between q and r for fiat universe models with mat- 
ter+quiessence is found by eliminating fix between equa- 
tion (14), with fi = 1, and equation (26). This gives 



r = 3{l + w)q- -(l + 3w), 



(27) 



which is the equation of the dotted straight lines in figure 
n When r^x = 1, all models lie on the solid curve given by 



3 1 

'^=r+2 

9 
r — -w{l + ■ 



or 



r = 2q^ 



(28) 
(29) 
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in accordance with equation (15) since q = for these 
models. This curve is the lower bound for all models with 
a constant w. For — 1 < w < 0, all matter -|-quiessence 
models will at any time fall in the sector between this 
curve and the r = 1-line which corresponds to ACDM. The 
results shown in Alam et al. (2003) seem to indicate that 
all matter -|-quintessence models will fall within this same 
sector as the matter -|-quiessence models do. However, as 
we will show below, this is not strictly correct. 

2.2. Scalar field models 

If the source of the dark energy is a scalar field (j>, as in 
the quintessence models (Wetterich 1988; Peebles & Ratra 
1988), the equation of state factor w is 



4'^ - 2^(0) 
02 + 2^(0) 

Differentiation gives 

2(i){24>V - <j)V) 
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2V 



Using the equation of motion of the scalar field 
dV 



(26) 0=-3iJ0 



(31) 



(32) 



(33) 



and V = (j)dV/ d(j) in equation (|31|) and inserting the result 
in equation (|24|l we obtain 



n + 127rG 



m 



SttG 



v_ 



(34) 



Evans et al.: Geometrical Constraints on Dark Energy 




Fig. 1. The q — r-plane for flat matter+quiessence mod- 
els. The horizontal curve has w = — 1 (ACDM). Then w 
increases by 1/10 counterclockwise until we reach w = 1 
in the upper right. When J7xo = all models start at the 
point q = 0.5, r = 1 (Einstein-de Sitter model). As fl^o in- 
creases every model moves towards the solid curve which 
marks n^a = 1- The crosses mark the present epoch. 



and furthermore, 

^ 3 «x 4^G /1.2 

Hence the statefinder s is 
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For models with matter -|-quintessence-|-curvature, the 
Friedmann and energy conservation equations give 
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As customary when discussing quintessence, we have in- 
troduced the Planck mass M^ — I/SttG. Furthermore, we 
have defined llkin = 0^/2pc, and f^pot = V{(j))/pc. For an 
exponential potential, V{(l)) — Aexp{—X(f>/M), and look- 
ing at values at the present epoch, one gets 
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2^m 
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Fig. 2. Present values of q and r for matter-|-quintessence 
with an exponential potential. 

Top panel: From top to bottom the different curves have 
A = -5, -4, -3, -2, -1, 0, 1, 2, 3, 4, 5. They aU start at the 
point qo(^kin — -73) — 1.595, r'o(rikin — -73) — 7.57) [mat- 
ter -fZeldovich gas (px = Px)]- As J^kin decreases when 
we move to the left, they join at the point (7o(^kin = 
0) = -0.595, ro(J7kin = 0) = 1) (ACDM, marked with 
a diamond). The dotted curve shows the area all mat- 
ter -fquiessence models must lie within at all times. 
Bottom panel: Zoom-in of the figure above. Here the curve 
having A^ = 2 is also plotted (thick line). 



Ehminating J7poto, using fimo + ^kino + ^ 
one obtains 



90 — -Z^ma — (1 — f^ko) + 3rikinO 



potO 



f^i-( 



fcO 



ro = (1 - r!ko) + 9^!: 



kinO 



3A-\/6r2kino(l — ^kO — f^mO — f^kino) 



(45) 
(46) 



By choosing for instance ^^q — 0.27 and f^ko = we 
can plot the values of q^ and tq for varying rikino! see 
figure [3 As we can see from equations (|45|l - (|46f) . when 
^kino — 0, go and tq are independent of A, and have the 
same values as in the ACDM model. This is obvious, since 
taking away the kinetic term will reduce quintessence to 
LIVE. However, when fikino is slightly greater then we 
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can make rg as large or as small as we like, by choosing 
|A| sufficiently large. There is no reason all quintessence 
models should lie inside the constant- w-curve. However, in 
order to get an accelerating universe today we must have 
A^ < 2. But also for A^ < 2 the present values of Qq and 
To can lie outside the constant-w-curve. In fact, when we 
move on to the time-evolving statefinders, plotting q and r 
as functions of time for given initial conditions, we obtain 
plots like figure|31 Here we have chosen as initial conditions 




-0.8 -0.6 -0.4 



0.2 0.4 



Fig. 3. Time-evolution of q and r for models with mat- 
ter and quintessence with an exponential potential. The 
crosses mark the present epoch. The diamond repre- 
sents the present ACDM model. The curve on top has 
A = 0.2 and then A increases by 0.2 for each curve go- 
ing counter-clockwise until we reach A = 2 to the right. 
The corresponding values for llkin today are fikino — 
0.002, 0.01, 0.02, 0.04, 0.06, 0.09, 0.12, 0.165, 0.22, 0.29. The 
dotted curve shows the area all matter-|-quiessence models 
must lie within at all times. We see that all models will 
eventually move towards this curve. 




Fig. 4. Present values of the statefinder parameters and 

the deceleration parameter for models with matter and 

quintessence with an exponential potential. The diamond 

represents the ACDM model. 

Top panel: From left to right the different curves have 

A = -5, -4, -3, -2, -1, 0, 1, 2, 3, 4, 5. 

Bottom panel: From top to bottom the different curves 

have A = -5, -4, -3, -2, -1, 0, 1, 2, 3, 4, 5. 



^mo = 0.27 and rj^o = as above, and h — 0.71. The last 
initial condition, for the quintessence field, we have chosen 
to be (f>o = A//100 combined with the overall constant A in 
the potential chosen to give rJkinO as stated in the caption 
of figure This corresponds to the universe being matter 
dominated at earlier times. When fipoto ^ f^kino we have 
high acceleration today. Choosing fikino = will again give 
us ACDM. The three rightmost curves in the figure have 
A^ > 2 and no eternal acceleration, although the A = 1.6 
universe accelerates today. It seems that in order to get a 
universe close to what we observe, r and q for models with 
matter -fquintessence with an exponential potential will 
essentially lie within the same area as matter-fquiessence 
models. In figure 01 we have plotted the trajectories in the 
SQ-ro-plane and the SQ-go-plane for the same models as in 
figure 121 to be compared with figures 5c and 5d in Alam 
et al. (2003). 



Choosing instead a power-law potential V{(j)) = Acf) " 
nves V — —^V and 
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(47) 
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We see that for (j)o = M we get the same curves in the 
(7o-ro-plane when varying a as we got when varying A 
in the exponential potential, see figure |2] We also see 
that varying ^q for a given value of a is essentially the 
same as varying a. Figure El shows the go^^o-plane for 
the case a = 2. Figure El shows an example of time- 
evolving statefinders {cpo = M, fikino = 0.05, Slmo = 0.27 
r^ko = 0, h=0.71). If one compares this plot with fig- 



ure lb in Alam et al. (2003) the two do not quite agree. 
Alam et al. (2003) do not give detailed information about 



the initial conditions for the quintessence field. Our initial 
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2.3. Dark energy fluid models 



Fig. 5. Present values of q and r for matter+quintessence 
with a power-law potential with a = 2. From 
top to bottom the different curves have 0o = 
8M,4M,2M,M, ^,^,^.The diamond represents the 
ACDM model. The dotted curve shows the area all mat- 
ter -|-quiessence models must lie within at all times. 



conditions correspond to a universe which was matter- 
dominated up to now, when quintessence is taking over. 




Fig. 6. Time-evolution of q and r for models with matter 
and quintessence with a power-law potential. The crosses 
mark the present epoch, the diamond represents the 
present ACDM model. All models start out from the hori- 
zontal ACDM line and will eventually end up as a deSitter 
Universe (g = —1,?" = 1). The curve going deepest down 
has a = 5 and moving upwards we have a = 4, 3, 2, 1. The 
dotted curve shows the area all matter-|-quiessence models 
must lie within at all times. Obviously the same is not the 
case for matter+quintessence models. 



We will now find expressions for r and s which are valid 
even if the dark energy does not have an equation of state 
of the formpx = wpx- This is the case e.g. in the Chaplygin 
gas models (Kamenshchik, Moschella & Pasquier 2001; 
Bilic, Tupper & VioUier 2002). The expression for the de- 
celeration parameter can be written as 
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3^]n. 



2 V Px 

and using this in equation H15|) we find 
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If the universe contains only dark energy with an equation 
of state p = p{p), then 
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. dp 



dp 

which leads to 
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(52) 

(53) 

(54) 



If the universe contains cold matter and dark energy these 
expressions are generalized to 



1 



9 Px + Px dp^ 
2 Pm + Px 9px 
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s = I 1 H 

Px 



dp^ 
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(55) 
(56) 



The Generalized Chaplygin Gas (GCG) has an equation 
of state (Bento, Bertolami & Sen 2002) 



A 
P" 



(57) 



and integration of the energy conservation equation gives 



P = 



A + Ba-3(i+") 



(58) 



where i? is a constant of integration. This can be rewritten 
as 



P = Po 



A, + (1 - A,)a:3(l+") 



(59) 



where po = {A+B)^/^^+'^\ and A^ ^ A/{A+B). For a flat 
universe with matter and a GCG, the Hubble parameter 
is given by 






a„ox=^+(i-a„o) ^ + (1 - A,)x3(i+") 



.(60) 
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This leads to the following expressions for q{x) and r{x): 2.5. The luminosity distance to third order in z 

•i O 3 _i_ n — O '\(^ — A \ f~^ The statefinder parameters appear when one expands 

q(x) — 5 — — — - — 57-5 1 (61) the luminosity distance to third order in the redshift z. 

2 ^moX'^ + (1 - ^mo) ■v'^l^ 



r(x) 
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■/(^) 



-/'(^), 



/i2(x)-" ' 2 h^{x) 

where (5 = 3(1 + a), h(x) — H{x)/Hq, and 

V = A, + {1- A,)x^ 
fix) ^ a„oa;2 + (l-0mo)(l-^)w' 



1 xP-\ 



(62) 



(63) 
(64) 



Although this expansion has been presented earlier (Chiba 
& Nakamura 1998; Visser 2003) in a slightly different no- 
tation, we carry out this derivation here for completeness. 
The luminosity distance is given by the expression 



dv = 



Ho^jW 



--Sk{^^\I), 



kO 



In the r-s plane, the GCG models will lie on curves given 
by (see Gorini, Kamenshchik & Moschella 2003) 



where 5k (x) = sin a; for fc = 1, S\^{x) — x for k 
5k (s) = sinhx for k — —1, and 



I = Hn 



^ 9s(s + l) 
2 a 



dz 



(65) 



We note that a recent comparison of GCG models with <^ / \ 
SNIa data found evidence for a > 1 (Bcrtolami et al. 
2004). 



H{^)' 
Using the approximation 
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one finds to third order in z 
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2.4. Cardassian models 
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As an alternative to adding a negative-pressure compo- Taylor expanding the Hubble parameter to second order 
nent to the energy-momentum tensor of the Universe, one in z, we have 
can take the view that the present phase of accelerated .z 1 

expansion is caused by gravity being modified, e.g. by the I Ri -, (74) 

presence of large extra dimensions. For a general discus- "^0 ^ / ^ 

sion of extra-dimensional models and statefinder parame- 
ters, see Alam & Sahni (2002). 

As an example, we will consider the Modified 
Polytropic Cardassian ansatz (MFC) (Freese & Lewis 
2002; Gondolo & Freese 2003), where the Hubble param- 
eter is given by 
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z=0 



H{x)=Ho\ln^ox^l 
with 



1/4, 



By using (1 + j/)" 



1 — y + y^ with y = jz + kz"^, we get 



(66) to third order in z 
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(77) 



and where n and ip are new parameters {tp is usually called 

q in the literature, but we use a different notation here to 

avoid confusion with the deceleration parameter). For this We wish to express 7 and k in terms of the deceleration 

model, the deceleration parameter is given by parameter 

"l + i 






1 



- 1 



(68) 



and the statefinder r by 



r{x) 



9 IJ 
4 1 



m(1 — n) — (1 + nu) 



2q 



(1 — n)^u 
(1 + m)(1 + nu)' 



(69) 



90 = 


1 aa\ 


and the statefinder 


ro = 


/a2a\ 

V d^ )t=to' 


From H = a/a one finds 


H = 


'{^ + q)H\ 


H = 


{r + iq + 2)H\ 



(78) 



(79) 



(80) 
(81) 
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From a = {1 + z) ^, one gets H = —z/{l + z), and hence 3.1. Model-independent constraints 



d 1 d 
dz {l + z)Hd 
d^ 1 d 


1 d' 


dz^ (1 + z)H dt 


{l + z)H dt 


After some algebra one then finds 


7 = 1 + 90, 

K = -(rn ~ ql). 





Substitution of these expressions for 7 and k in 



(82) 
(83) 

(84) 
(85) 

l(77|) gives 



The approximation to di^ in equation (|87|l is independent 
of cosmological model, the only assumption made is that 
the Universe is described by the Friedmann-Robertson- 
Walker metric. We see that, in addition to Hq, this third- 
order expansion of c^l depends on qq and the combination 
rp — fiko- Fitting these parameters to the data, we find the 
constraints shown in Fig.Q The results are consistent with 
those of similar analyses in Caldwell & Kamionkowski 
(2004) and Gong (2004). In figure|Hlwe show the marginal- 



(1 + z)/ w z 



l+2(l-'?o)^- g(l+ro 



qo - 3qo)z^ 



and equation H73|l then finally leads to 



(86) 



a 2- 



dL « — 
-no 



1. N 1/ 




3(7o - ^ko)z'^ 



(87) 



Fig. 7. Likelihood contours (68, 95 and 



resulting 



r\ 1 c J ■ c j-i. i 1 from a fit of the expansion of the luminosity distance to 

Une can also find an expression tor the present value . . ^ ■' 

of the time derivative of the equation of state parameter 

w in terms of the statefinder rg. A Taylor expansion to 

first order in z gives 

w{z) = Wo 7]—z, 

From equation (|14|l we have 

_ 2go — ^mO ~ ^xO 
3f^xo 

and from equation (|24|l we get 



(89) 



^(*o) = o7T^"i ^mO 
6 StxO 



1 + ^^"o(l + Wo) 



nxo-roL(90) 



Hence 



w{z) wwo - - 



,9 f^niO - To 

1 + :^wo(l + Wo) + — -^ 



-113 5 



(91) 



3. Lessons drawn from current SNIa data 



Fig. 8. Marginalized probability distributions for go (full 
line) and r^ — J7ko (dotted line). 



In this section we will consider the SNIa data presently 
available, in particular whether one can use them to learn 
about the statefinder parameters. We will use the recent 
collection of SNIa data in Riess et al. (2004), their 'gold' 
sample consisting of 157 supernovae at redshifts between 
~ 0.01 and - 1.7. 



ized distributions for qo and ro — r^ko- We note that the 
supernova data firmly support an accelerating universe, 
go < at more than 99% confidence. However, about 
the statefinder parameter r^, little can be learned with- 
out an external constraint on the curvature. Imposing a 
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flat universe, e.g. by inflationary prejudice or by invok- 
ing the CMB peak positions, there is still a wide range 
of allowed values for rg. This is an indication of the lim- 
ited ability of the current SNIa data to place constraints 
on models of dark energy. There is only limited informa- 
tion on anything beyond the present value of the second 
derivative of the Hubble parameter. 

Under the assumption of a spatially flat universe, 
^ko = 0, with flmo — 0.3, one can use equations H89|l 
and (|91|l to obtain constraints on wq and wi in the ex- 
pansion w{z) = Wq + wiz of the equation of state of dark 
energy. The resulting likelihood contours are shown in fig- 
ure|51 As can be seen in this figure, there is no evidence for 
time evolution in the equation of state, the observations 
are consistent with wi = 0. The present supernova data 
show a slight preference for a dark energy component of 
the 'phantom' type with wq < —I (Caldwell 2002). Note, 
however, that the relatively tight contours obtained here 
are caused by the strong prior Jlmo = 0.3. It should also be 
noted that the third-order expansion of di, is not a good 
approximation to the exact expression for high z and in 
some regions of the parameter space. 



Table 1. Results of fitting the models considered in this 
subsection to the SNIa data 




Fig. 9. Likelihood contours (68, 95 and 99%) for the co- 
efficients Wq and wi in the linear approximation to the 
equation of state w(z) of dark energy, resulting from a fit 
of the expansion of the luminosity distance to third order 



3.2. Direct test of models against data 

The standard way of testing dark energy models against 
data is by maximum likelihood fitting of their parameters. 
In this subsection we will consider the following models: 

1. The expansion of d]^ to second order in z, with h and 
go as parameters. 

2. The third-order expansion of d]^, with h, qo, and rg — 
r^ko as parameters. 

3. Flat ACDM models, with fimo and h as parameters to 
be varied in the fit. 



Model 


Anim 


^ parameters 


B 


2. order expansion 


177.1 


2 


187.2 


3. order expansion 


162.3 


3 


177.5 


Flat ACDM 


163.8 


2 


173.9 


ACDM with curvature 


161.2 


3 


176.4 


flat -1- constant EoS 


160.0 


3 


175.2 


MPC 


160.3 


4 


180.5 


GCG 


161.4 


4 


181.6 


Alam et al. 


160.5 


4 


180.7 



ACDM with curvature, so that ilmO: ^Ao (the contri- 
bution of the cosmological constant to the energy den- 
sity in units of the critical density, evaluated at the 
present epoch), and h are varied in the fits. 
Flat quiessence models, that is, models with a constant 
equation of state w for the dark energy component. 
The parameters to be varied in the fit are ri,nO, w, and 
h. 

The Modified Polytropic Cardassian (MPC) ansatz, 
with f2mo, q, n, and h as parameters to be varied. 
The Generalized Chaplygin Gas (GCG), with rimOi ^s, 
a, and h as parameters to be varied. 
The ansatz of Alam et al. (2003), 



H = HoVn„ 



iOX" 



Ao + Aix + A2x'^, 



(92) 



where we restrict ourselves to flat models, so that Aq = 
1 — rimo — Al — A2. The parameters to be varied are 
rimo, Al, A2, and h. 

Note that these models have different numbers of free 
parameters. To get an idea of which of these models 
is actually preferred by the data, we therefore employ 
the Bayesian Information Criterion (BIG) (Schwarz 1978; 
Liddle 2004). This is an approximation to the Bayes fac- 
tor (Jeffreys 1961), which gives the posterior probability 
of one model relative to another assuming that there is 
no objective reason to prefer one of the models prior to 
fltting the data. It is given by 



■8==Xmi„ + ^parln/^data, 



(93) 



where Xmin i^ the minimum value of the x^ for the given 
model against the data, A^par is the number of free param- 
eters, and A^data is the number of data points used in the 
flt. As a result of the approximations made in deriving it, 
B is given in terms of the minimum x^, even though it is 
related to the integrated likelihood. The preferred model 
is the one which minimizes B. In tableHwe have collected 
the results for the best-fltting models. When comparing 
models using the BIG, the rule of thumb is that a differ- 
ence of 2 in the BIG is positive evidence against the model 
with the larger value, whereas if the difference is 6 or more, 
the evidence against the model with the larger BIG is con- 
sidered strong. The second-order expansion of c^l is then 
clearly disfavoured, thus the current supernova data give 
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information, although limited, on tq — ^ko- We see that 
there is no indication in the data that curvature should be 
added to the ACDM model. Also, the last three models 
in tabled seem to be disfavoured. One can conclude that 
there is no evidence in the current data that anything be- 
yond flat ACDM is required. This does not, of course, rule 
out any of the models, but tells us that the current data 
are not good enough to reveal physics beyond spatially 
flat ACDM. A similar conclusion was reached by Liddle 
(2004) using a more extensive collection of cosmological 
data sets and considering adding parameters to the flat 
ACDM model with scale- invariant adiabatic fluctuations. 

3.3. Statefinder parameters from current data 

If the luminosity distance ^l is found as a function of 
redshift from observations of standard candles, one can 
obtain the Hubble parameter formally from 



Hix) 



d_ 
dx 



n -1 



(94) 



However, since observations always contain noise, this re- 
lation cannot be applied straightforwardly to the data. 
Alam et al. (2003) suggested parametrizing the dark en- 
ergy density as a second-order polynomial in x, px = 
yOco(^o + Aix + A2X^), leading to a Hubble parameter of 
the form 



H{x) = Ho^/n^ox^ +Ao + Aix + A2X^, 



(95) 



and fitting Aq, Ai, and A2 to data. This parametrization 
reproduces exactly the cases w = —1 {Ai = A2 =0), 
w = -2/3 (Ao =A2= 0), and w = -1/3 [Aq ^ Ai ^ 0), 
and the luminosity distance-redshift relationship is given 

by 



l + z 
Ho 



1+2 



dx 



vf^moa 



Ao + Aix - 



(96) 



Having fitted the parameters Ao, Ai, and A2 to e.g. su- 
pernova data using H96|l . one can then find q and r by 
substituting equation (pl5|l into (0) and ©: 



q{x) 
r{x) 



1- 



2Aix + 3Ao 






- A2X^ 

Ao 



Aix + Ao 



n„,oX^ + A0+A1X + ^2X2 ' 

and furthermore the statefinder s is found to be 

, - _ 2 Aix + A2X^ 
^^"^^ ~ 3 3Ao + 2Aix + A2X^ ' 

and the equation of state is given by 

1 Aix + 2A2x^ 
'w(x) = —1 H 



3 Ao + Aix + A23 



(97) 
(98) 

(99) 
(100) 



The simulations of Alam ct al. (2003) indicated that the 
statefinder parameters can be reconstructed well from sim- 
ulated data based on a range of dark energy models, 



so we will for now proceed on the assumption that the 
parametrization in equation (|95() is adequate for the pur- 
poses of extracting q, r and s from SNIa data. We com- 
ment this issue in section 4. 

In figure 1101 we show the deceleration parameter q 
and the statefinder r extracted from the current SNIa 
data. The error bars in the figure are la limits. We have 
also plotted the model predictions for the same quantities 
(based on best-fitting parameters with errors) for ACDM, 
quiessence, and the MPC. The figure shows that all models 
are consistent at the la level with q and r extracted us- 
ing equation (|95|) . Thus, with the present quality of SNIa 
data, the statefinder parameters are, not surprisingly, no 
better at distinguishing between the models than a direct 
comparison with the SNIa data. We next look at simu- 
lated data to get an idea of how the situation will improve 
with future data sets. 

4. Future data sets 

We will now make an investigation of what an idealized 
SNIa survey can teach us about statefinder parameters 
and dark energy, following the procedure in Saini, Weller 
& Bridle (2004). 

A SNAP-like satellite is expected to observe ^2000 
SN up to z '^ 1.7. Dividing the interval < a; < 1.7 
into 50 bins, we therefore expect ~40 observations of SN 
in each bin. Empirically, SNIa are very good standard 
candles with a small dispersion in apparent magnitude 
Cmag = 0.15, and there is no indication of redshift evolu- 
tion. The apparent magnitude is related to the luminosity 
distance through 



m{z) =M + 5logDi,{z), 



(101) 



where M = Mq + 51og[iJo ^Mpc ^] + 25. The quantity 
Mo is the absolute magnitude of type la supernovae, and 
Dhiz) — Hodi^{z) is the Hubble constant free luminosity 
distance. The combination of absolute magnitude and the 
Hubble constant, Ai, can be calibrated by low-redshift 
supernovae (Hamuy et al. 1993; Perlmutter et al. 1999). 
The dispersion in the magnitude, CTmag, is related to the 
uncertainty in the distance, a, by 



In 10 



dhiz) 



(102) 



and for CTmag = 0.15, the relative error in the luminosity 
distance is ~ 7%. If we assume that the d^ we calculate 
for each z value is the mean of all d^s in that particular 
bin, the errors reduce from 7% to 0.07/\/40 w 0.01 = 1%. 
We do not add noise to the simulated dj^, and hence our 
results give the ensemble average of the parameters we fit 
to the simulated data sets. 

4.1. A ACDM universe 

We first simulate data based on a flat ACDM model with 
^mo = 0.3, h = 0.7, giving the data points shown in figure 
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Fig. 10. The deceleration parameter q and the statefinder 
r extracted from current SNIa data using the Alam 
parametrization of H (top row), for ACDM (second 
row) , quiessence (third row) , and the Modified Polytropic 
Cardassian ansatz (bottom row) 



111! To this data set we first fit the quiessence model, the 
MPC, the GCG, and the parametrization of H from equa- 
tion H95() . Since aU models reduce to ACDM for an appro- 
priate choice of parameters, distinguishing between them 
based on the x^ per degree of freedom alone is hard. Based 
on the best-fitting values and error bars on the parame- 
ters Aq, Al, and A2 in equation H95|) we can reconstruct 
the statefinder parameters from eqs. H97I) - (|99|l . In figure 
[T^ we show the deceleration parameter and statefinder 
parameters reconstructed from the simulated data. The 
statefinders can be reconstructed quite well in this case, 



Fig. 11. Binned, simulated data set for a flat ACDM uni- 
verse with rimo = 0.3. The la error bars are also shown. 



1.0 r 





0,0 0.5 1.0 1.5 2.0 



0,0 0.5 i.i 



-0.6 -0.4 -0,2 



Fig. 12. The statefinder parameters and the deceleration 
parameter for the best-fitting reconstruction of the simu- 
lated data based on ACDM, using the parametrization of 
Alam et al. The Ict error bars are also shown. 



e.g. we see clearly that r is equal to one, as it should for 
flat ACDM. In figure ^| we show the statefinders for a 
selection of models, obtained by fitting their respective 
parameters to the data, and using the expressions for q 
and r for the respective models derived in earlier sections, 
e.g. equation (|61l) and (|62() for the Chaplygin gas. Since all 
models reduce to ACDM for the best-fitting parameters, 
their q and r values are also consistent with ACDM. Thus, 
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Fig. 13. The statefinder parameters for a selection of Fig. 14. Simulated luminosity distance-redshift relation- 
models, evaluated at the best-fitting values of their re- ship for the Generalized Chaplygin Gas with A^ = 0.4, 
spective parameters to the simulated ACDM datasct, with a = 0.7. 
errors included. 



if the dark energy really is LIVE, a SNAP-type experiment 
should be able to demonstrate this. 



4.2. A Chaplygin gas universe 

We have also carried out the same reconstruction exercise 
using simulated data based on the GCG with As = 0.4, 
a = 0.7, see figure [T^ Figure [T31 shows q and r recon- 
structed using the parametrization of H. The same quan- 
tities for the models considered, based on their best-fitting 
parameters to the simulated data, are shown in figure El 
For the Cardasslan model, the best-fitting value for the 
parameter n, ribf, depends on the extent of the region 
over which we allow n to vary. Extending this region to 
larger negative values for n moves nbf in the same direc- 
tion. However, the minimum x^ value does not change 
significantly. This is understandable, since H{x) for the 
MFC model is insensitive to n for large, negative values 
of n. The quantities r{x) and q{x) also depend only weakly 
on the allowed range for n, whereas their error bars are 
sensitive to this parameter. We chose to impose a prior 
n > — 1, producing the results shown in figure ITHl The 
best-fitting values for tp and n were, respectively, 0.06 and 
-0.94. 

Figure [T7I shows the deceleration parameter extracted 
from the Alam et al. parametrization (full line), with 
Icr error bars. Also plotted is the best fit q{z) from 
the quiessence (squares), Cardasslan (triangles) and 
Chaplygin (asterisk) models. We note that the q{z) from 
the Alam et al. parametrization has a somewhat deviat- 
ing behaviour from the input model, especially at larger 
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Fig. 15. The statefinder parameters and the deceleration 
parameter for the best-fitting reconstruction of the simu- 
lated data based on the GCG, using the parametrization 
of Alam et al. The Itr error bars are also shown. 



z. Also, no model can be excluded on the basis of their 
predictions for q{z) 

Figure [THl shows the same situation for the statefinder 
parameter r(z). Note again that for large z, the recovered 
statefinder from the Alam et al. parametrization does not 
correspond well with the input model. As with the case 
for q{z), the quiessence and Cardasslan models follow each 
other closely. These, however, do not agree with the input 
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Fig. 16. The statefinder parameters for a selection of 
models, evaluated at the best-fitting values of their re- 
spective parameters to the simulated Chaplygin gas data 
set, with errors included. 



Fig. 18. Comparison of r{z) extracted using the 
parametrized H(z) with r{z) for the various best-fitting 
models. Error bars are only shown on the values extracted 
using the Alam et al. parametrization, but in the other 
cases they are roughly of the same size as the symbolsS. 



model for low values of z (similar to the case for q{z) they 
diverge for low z). Comparing the statefinder r for the 
quiessence and Cardasslan models with that of the input 
GCG model, indicates that, not surprisingly, neither of 
them is a good fit to the data. 
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Fig. 17. Comparison of q{z) extracted using the 
parametrized H{z) with q{z) for the various best-fitting 
models. Error bars are only shown on the values extracted 
using the Alam et al. parametrization, but in the other 
cases they are roughly of the same size as the symbols. 



This exercise indicates that the statefinders can po- 
tentially distinguish between dark energy models, if q, r 
and s can be extracted from the data in a reliable, model- 
independent way. However, the fact that r extracted from 
the simulated data using the Alam et al. parametrization 
docs not agree well with the true r of the underlying model 
for z > 0.5, indicates that one needs a better parametriza- 
tion in order to use statefinder parameters as empirical 
discriminators between dark energy models. 

5. Conclusions 

We have investigated the statefinder parameters as a 
means of comparing dark energy models. As a theoreti- 
cal tool, they are very useful for visualizing the behaviour 
of different dark energy models. Provided they can be ex- 
tracted from the data in a reliable, model-independent 
way, they can give a first insight into the type of model 
which is likely to describe the data. However, SNIa data of 
quality far superior to those presently available are needed 
in order to distinguish between the different models. And 
even with SNAP-quality data, there may be difficulties in 
distinguishing between models based on the statefinder 
parameters alone. Furthermore, the parametrization of 
H{z) used here and in Alam et al. (2003) is probably not 
optimal, as shown in section 4.2. The same conclusion was 
reached in a recent investigation by Jonsson et al. (2004), 
where they considered reconstruction of the equation of 
state w{x) from SNfa data using equation (|100|l . They 
found that this parametrization forces the behaviour of 
w{x) onto a specific set of tracks, and may thus give spuri- 
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ous evidence for redshift evolution of the equation of state. 
Although this conclusion has been contested by Alam et 
al. (2004), it is clear that finding a parametrization which 
is sufficiently general, and at the same time with reason- 
ably few parameters is an important task for future work. 
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